The high mutation rate at the D614G hotspot-furin cleavage site region increases the priming efficiency of the Spike protein by furin protease: analysis of Indonesian SARS-CoV-2 G614 variants obtained during the early COVID-19 pandemic.

D614G mutation plays a significant role in the transmissibility of SARS-CoV-2. Identification of other mutations related to D614G mutation within the Spike protein is pivotal as they might contribute to the pathogenicity of SARS-CoV-2. This study aims to analyze the mutation rate of furin cleavage site (FCS) region of Indonesian origin SARS-CoV-2 and to predict the effect of mutation against Spike priming efficiency by furin. A total of 375 sequences of Indonesian isolates obtained during the early pandemic were used for mutation analysis. Mutation analysis includes mutation pattern, variability, frequency of mutation, amino acid conservation, and mutation rate. The effect of mutation against Spike priming efficiency by furin protease from eight sequences with mutation in the FCS region was analyzed by protein-protein docking. We showed that mutations related to the G614 variant were increasing through time, in contrast to the D614 variant. The FCS region at the position 675-692 contained the most variable (66.67%) as well as the highest mutation frequency (85.92%) and has been observed to be the hotspot mutations linked to the D614G mutation. The D614G hotspot-FCS region (residue 600-700) had the highest amino acid change per site (20.8%) as well as the highest mutation rate as 1.34 × 10-2 substitution per site per year (95% CI 1.79 × 10-3-2.74 × 10-2), compared with other Spike protein regions. Mutations in the FCS region were the most common mutation found after the D614G mutation. These mutations were predicted to increase the Spike priming efficiency by furin. Thus, this study elucidates the importance of D614G mutation to other mutations located in the FCS region and their significance to Spike priming efficiency by furin. Supplementary Information: The online version contains supplementary material available at 10.1007/s13337-023-00827-w.

Co-infection of SARS-CoV-2 with other viral respiratory pathogens in Yogyakarta, Indonesia: A cross-sectional study.

Background: Growing evidence shows that viral co-infection is found repeatedly in patients with Coronavirus Disease-2019 (COVID-19). This is the first report of SARS-CoV-2 co-infection with viral respiratory pathogens in Indonesia. Methods: Over a one month period of April to May 2020, SARS-CoV-2 positive nasopharyngeal swabs in our COVID-19 referral laboratory in Yogyakarta, Indonesia, were tested for viral respiratory pathogens by real-time, reverse transcription polymerase chain reaction (RT-PCR). Proportion of co-infection reported in percentage. Results: Fifty-nine samples were positive for other viral respiratory pathogens among a total of 125 samples. Influenza A virus was detected in 32 samples, Influenza B in 16 samples, Human metapneumovirus in 1 sample, and adenovirus in 10 samples. We did not detect any co-infection with respiratory syncytial virus. Nine (7.2%) patients had co-infection with more than two viruses. Conclusion: Viral co-infection with SARS-CoV-2 is common. These results will provide a helpful reference for diagnosis and clinical treatment of patients with COVID-19.

Evolutionary dynamics of SARS-CoV-2 circulating in Yogyakarta and Central Java, Indonesia: sequence analysis covering furin cleavage site (FCS) region of the spike protein.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a new virus responsible for the COVID-19 pandemic. The emergence of the new SARS-CoV-2 has been attributed to the possibility of evolutionary dynamics in the furin cleavage site (FCS) region. This study aimed to analyze the sequence of the FCS region in the spike protein of SARS-CoV-2 isolates that circulated in the Special Region of Yogyakarta and Central Java provinces in Indonesia. The RNA solution extracted from nasopharyngeal swab samples of confirmed COVID-19 patients were used and subjected to cDNA synthesis, PCR amplification, sequencing, and analysis of the FCS region. The sequence data from GISAID were also retrieved for further genome analysis. This study included 52 FCS region sequences. Several mutations were identified in the FCS region, i.e., D614G, Q675H, Q677H, S680P, and silent mutation in 235.57 C > T. The most important mutation in the FCS region is D614G. This finding indicated the G614 variant was circulating from May 2020 in those two provinces. Eventually, the G614 variant totally replaced the D614 variant from September 2020. All Indonesian SARS-CoV-2 isolates during this study and those deposited in GISAID showed the formation of five clade clusters from the FCS region, in which the D614 variant is in one specific cluster, and the G614 variant is dispersed into four clusters. The data indicated there is evolutionary advantage of the D614G mutation in the FCS region of the spike protein of SARS-CoV-2 circulating in the Special Region of Yogyakarta and Central Java provinces in Indonesia.

Detection of SARS-CoV-2 spike protein D614G mutation by qPCR-HRM analysis.

OBJECTIVES: Monitoring the spread of the G614 in specific locations is critical as this variant is highly transmissible and can trigger the emergence of other mutations. Therefore, a rapid and accurate method that can reliably detect the D614G mutation will be beneficial. This study aims to analyze the potential use of the two-step Reverse Transcriptase quantitative polymerase chain reaction - high resolution melting analysis (RT-qPCR-HRM) to detect a specific mutation in the SARS-CoV-2 genome. METHODS: Six SARS-CoV-2 RNA samples were synthesized into cDNA and analyzed with the qPCR-HRM method in order to detect the D614G mutation in Spike protein of SARS-CoV-2. The primers are designed to target the specific Spike region containing the D614G mutation. The qPCR-HRM analysis was conducted simultaneously, and the identification of the SARS-CoV-2 variant was confirmed by conventional PCR and Sanger sequencing methods. RESULTS: The results showed that the melting temperature (Tm) of the D614 variant was 79.39 ± 0.03 °C, which was slightly lower than the Tm of the G614 variant (79.62 ± 0.015 °C). The results of the HRM analysis, visualized by the normalized melting curve and the difference curve were able to discriminate the D614 and G614 variant samples. All samples were identified as G614 variants by qPCR-HRM assay, which was subsequently confirmed by Sanger sequencing. CONCLUSIONS: This study demonstrated a sensitive method that can identify the D614G mutation by a simple two-step RT-qPCR-HRM assay procedure analysis, which can be useful for active surveillance of the transmission of a specific mutation.